FROM LIMIT CYCLES TO STRANGE ATTRACTORS 
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Abstract. We define a quantitative notion of shear for limit cycles of flows. We prove that strange attrac- 
tors and SRB measures emerge when systems exhibiting limit cycles with sufficient shear are subjected to 
periodic pulsatile drives. The strange attractors possess a number of precisely-defined dynamical properties 
that together imply chaos that is both sustained in time and physically observable. 
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1. Introduction 

This paper is about a mechanism for producing chaos: shear. We are guided by the idea that in the 
presence of shear, a stable dynamical structure can be transformed into a strange attractor with strong 
stochastic properties by forcing the structure with a pulsatile drive. The forcing does not overwhelm the 
intrinsic dynamics. Instead, it acts as an amplifier, amplifying the effects of the intrinsic shear. We focus 
on one particular dynamical structure of great importance: the limit cycle. Limit cycles are asymptotically 
stable periodic orbits of flows on Riemannian manifolds. 

The application of a periodic pulsatile drive to a flow exhibiting a limit cycle causes deformations to occur. 
If shear is present in a neighborhood of the limit cycle, if the limit cycle only weakly attracts nearby orbits, 
and if the time between pulses (the relaxation time) is sufficiently large, then stretch-and-fold geometry 
emerges in a neighborhood of the limit cycle. Stretch-and-fold geometry suggests that chaotic behavior that 
is both sustained in time and observable may exist. We prove that such chaotic behavior does exist in a 
certain parameter regime for any (generic) forcing function if the shear is sufficiently strong. Moreover, we 
define a quantity called the shear integral that quantifies the amount of shear that is present in the intrinsic 
flow in a neighborhood of the limit cycle. We emphasize that the shear integral depends only on the intrinsic 
system and not on the external forcing. Our result is the first of its kind for general limit cycles. Wang and 
Young [16, 17] obtain results of a similar flavor for supercritical Hopf bifurcations and certain linear models. 

The search for and analysis of stochastic behavior in deterministic dynamical systems have played a major 
role in guiding dynamical systems research. We discuss a few relevant developments. The theory of uniformly 
hyperbolic systems is well-developed. Let M be a compact Riemannian manifold and let / : M — > M be a 
C 2 diffeomorphism of M. An attractor for / is a compact set fi satisfying /(fi) = fi for which there exists 
an open set U C M (the basin) such that f(U) C U and £1 = Hi^o An attractor f2 is said to be an 

Axiom A attractor if the tangent bundle over splits into 2 D /-invariant subbundles E 8 and E u such that 
vectors in E s are contracted by Df and vectors in E u are expanded by Df (we assume E u is nontrivial) . An 
Axiom A attractor supports a special invariant measure known as a Sinai-Ruelle-Bowen (SRB) measure that 
describes the asymptotic distribution of the orbit of almost every point in U with respect to Riemannian 
volume and has strong stochastic properties. In this sense, the chaotic behavior associated with Axiom A 
systems is observable. It is also sustained in time because of the presence of positive Lyapunov exponent (s). 
One can, in principle, detect the presence of uniform hyperbolicity in a given system by finding invariant 
cone families with suitable properties. For example, Tucker uses this approach to prove that the Lorenz 
equations are chaotic for the classical parameter values studied by Lorenz [13]. 

Many systems of interest in the biological and physical sciences display some form of hyperbolicity but are 
not uniformly hyperbolic. A mature theory of nonuniform hyperbolicity has emerged over the last 4 decades. 
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However, the following problem remains a challenge. Given a dynamical system (or a parametrized family 
of dynamical systems), how can nonuniform hyperbolicity be detected? Numerical techniques include the 
calculation of Lyapunov exponents and the 0-1 test [4, 5]. This paper addresses the analytical component 
of the problem in the context of limit cycles. Our proofs are based on the recently-developed theory of rank 
one maps [15, 18]. Rank one theory is based on the ideas of Jakobson [8], Benedicks and Carleson [1, 2], and 
Young [19, 20]. Rank one theory provides checkable conditions that imply the existence of SRB measures 
with strong stochastic properties in parametrized families of diffeomorphisms. 

We conclude the introduction with a remark that the results obtained in this paper are in some sense dual 
to the phenomenon known as self-induced stochastic resonance (SISR) (see e.g. [3]). Our results demonstrate 
that certain intrinsic characteristics of a deterministic system (shear) can produce stochastic-type behavior 
when the system is forced in a deterministic way. SISR demonstrates that underlying phase space structures 
can produce deterministic (coherent) behavior in stochastically-forced systems when the noise level is taken 
to along certain distinguished limits. 

2. Statement of results 

We state the main results and discuss their relationship to the existing literature. Let / : E™ — > 1" be a 
C 5 vector field and consider the differential equation 

(2.D f -/(.). 

We assume that (2.1) admits an asymptotically stable hyperbolic periodic solution r] of length L and period 
Pq. Let 7 : M — > K™ be a function of the parameter s that parametrizes r\ by length. Define L = {-y(s) : 
s G [0,L)}. Solutions to (2.1) that begin sufficiently close to T will converge to T at an exponential rate as 
t — > oo. We are interested in the effects of adding periodic pulsatile forcing to the vector field defining (2.1). 
For < p < T, define the periodic function P p ,t : K — >■ K as follows. For ^ t ^ T, set 



1, HO^t^p 
0, if p < t < T 



and then extend periodically to all i € R by requiring P p .x(i + T) = P p T(t). We study the externally- forced 
system 

(2-2) —=f(x)+eP p , T (t)F(x) 

where F : K" — > R™ is a C 4 vector field and the parameter e > controls the amplitude of the forcing. 
Notice that the right side of (2.2) is not continuous. 

In Section 3 we compute a normal form of equation (2.2) that is valid in a tubular neighborhood M f=a VxD, 
where D is a closed disk in M™ -1 of sufficiently small radius. We are interested in the dynamics of (2.2) in 
the tubular neig hborhood M « T x\D. Since the external forcing is periodic with period T, it is natural 
to study the time-T map induced by (2.2). We write the time-T map as the composition of a kick map 
Hk : M — > M and a relaxation map H r : M — > int(M). Let Hk be the time-p map induced by the flow 
associated with (2.2). Notice that the external forcing is active during the kick phase because P p ^r{t) = 1 
for ^ t ^ p. For e sufficiently small, Hk maps M into M diffeomorphically. Let H r be the time-(T — p) 
map induced by (2.2) with e set to 0. There exists To = Tq(e) such that if T ^ To, then H r maps M into 
int(M). The composition Gt '■ = H r o Hk is the time-T map induced by (2.2). 

The dynamical properties of Gt ■ M — > int(M) depend on a number of factors. One feature common to 
every map Gt for T T is the existence of an attractor f2 defined by 

oo 

0= f]G T (M). 

i=0 

We call U := int(M) the basin of attraction of fi. For every x € U, G T (x) —> 57 as i — > oo. Two 
characteristics of the intrinsic system (2.1) play a key role in determining the structure of and the dynamical 
properties of Gt- shear and the strength of the limit cycle. We quantify these notions momentarily; for 
now, imagine that (2.1) exhibits strong shear in M if for most points x £ T, the velocity vector f(x) varies 
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substantially as x moves away from x in directions orthogonal to the limit cycle T. Think of the limit cycle 
r as strongly stable if solutions to (2.1) that begin in M converge quickly to V. If the shear is weak and the 
limit cycle is strongly stable, then the attractor fi associated with Gt will be an invariant closed curve. We 
are interested in the opposite situation. Suppose that the shear is strong in M and the limit cycle is weakly 
stable. The addition of the periodic pulsatile external force eP Pi T(t)F(x) will amplify the effect of the shear 
in the following way: disturbances that are created when P p ,t = 1 will be stretched during the relaxation 
period (when P Pt T = 0). The stretching effect increases in intensity as T increases. If T is large, then folds 
will be created in the phase space. If Gt exhibits stretch-and-fold geometry, then Gt potentially exhibits 
chaotic behavior that is sustained in time and observable. 
This paper aims to accomplish the following. 

(1) We define a computable quantity called the shear integral that quantifies the shear associated with the 
intrinsic system (2.1) near the limit cycle T. 

(2) We prove that if the magnitude of the shear integral is sufficiently large and if the contraction near the 
limit cycle T is sufficiently weak, then the following holds for suitable values of e. For a typical external 
vector field F, there exists 1\ > and a set A C [T^oo) of positive Lebesgue measure such that for 
T € A, the time-T map Gt associated with (2.2) admits a strange attractor fl and exhibits chaos that 
is sustained in time and observable. 

The quantity T\ satisfies T\ 3> p, ensuring sufficient relaxation time for the stretch-and-fold geometry to 
emerge. The term strange attractor refers to a number of precisely defined dynamical and structural 
properties that represent sustained, observable chaos. For T € A, Q supports a unique ergodic SRB mea- 
sure v. Here the term SRB measure refers to a measure v with a positive Lyapunov exponent v almost 
everywhere and whose conditional measures on unstable manifolds are absolutely continuous with respect to 
Riemannian volume on these manifolds. The SRB measure v satisfies the central limit theorem and exhibits 
exponential decay of correlations for Holder continuous observables. For Lebesgue almost every x in the 
basin of attraction U, the orbit of x has a positive Lyapunov exponent and is asymptotically distributed 
according to v in the sense that for every continuous function ip : U — > K, we have 

TO— 1 „ 

(2.3) lim — S^<p{G T {x))= ydv. 

TO— fOO JJT ' / 

i=0 J 

Notice that this statement is substantially stronger than the conclusion of the Birkhoff ergodic theorem. The 
Birkhoff ergodic theorem implies that (2.3) holds for v almost every x. However, v is singular with respect 
to Lebesgue measure (supported on a set of Lebesgue measure zero) because the dynamics are dissipative. 
We prove that (2.3) holds for Lebesgue almost every x E U. See (SA1)~(SA4) in Section 4 for a more 
precise description of the dynamical properties of Gt for T € A. 

We now define the shear integral. In Section 3 we derive a normal form of (2.1) that is valid in M . The 
normal form, expressed in the natural (s, ,z)-coordinates introduced in Section 3.1, is given by 

(2.4a) ~ = ||/( 7 ( 5 ))ir 1 + (P(s), z) + ^(s, z) 

as 

(2.4b) ~ = Az +u> 2 (s,z) 

as 

Here (•, •) denotes the inner product on Functions depending on s in (2.4a)-(2.4b) are periodic in s 

with period 2L. The matrix A is in Jordan canonical form. The functions wi and u>2 represent higher order 
corrections. The function f3 gives the pointwise magnitude and direction of the shear. Define the shear 
integral S by 

E = (Ei,...,S n _i) := / (3{r)dT 
Jo 

and define the shear factor a by a := ||S||. 

Having defined the shear integral, we describe the setting of the main theorem. We identify intrinsic 
parameters (parameters associated with /) and external parameters (parameters associated with the external 
forcing). We fix the normalized shear vector — and view the shear factor a as the first intrinsic parameter. 
The second intrinsic parameter quantifies the strength of the contraction near the limit cycle and is derived 
from A. We assume for the sake of simplicity that A is a diagonal matrix given by A = diag(Ai, . . . , A n _i) 

3 



^ \ n—1 
L'iPi 



where > Ai ^ A 2 > • • • ^ A„_i are the eigenvalues of A. We fix the eigenvalue ratios pn = j 1 for 
1 ^ i ^ n — 1 and we view the weakest eigenvalue Xi as an intrinsic parameter. The only external parameter 
is e, the factor that controls the amplitude of the external forcing. We fix p > 0. A key parameter derived 
from e, a, and Ai is the hyperbolicity factor t^j. 

One additional ingredient is needed. Even if a is large and |Ai| is small, a strange attractor cannot emerge 
unless the forcing F acts in direction(s) in which shear is present. We express this idea by introducing a 
certain function on the circle S := -^g- We identify S with the interval [0, 2L). In Section 3 we derive a 
normal form of the forced system (2.2) that is valid in M when the forcing is active (P p ,t = 1) : 

(2.5a) §- = ||/( 7 ( 5 ))ir 1 + (P(s), z) + u> s (s, z) 

as 

(2.5b) ^ = Az +eC(s) + uJs,z) 

as 

Functions depending on s in (2.5a)-(2.5b) are periodic in s of period 2L. The functions W3 and W4 are higher 
order corrections. The function £ is related to the projection of F in directions orthogonal to T. For s £ S, 
define s implicitly by 

P= f\\f(l(r))\\- l dr. 

J s 

Define the vector 

d : = 

and define $ : § — !> K by 

(2.6) $(s Q ) = (d, j\(r)dr\. 

We say that $ is a Morse function if the critical set C($) = {s £ S : $'(s) = 0} is finite and if for every 
s € C($>), we have <£>"(s) 7^ 0. We are now in position to state the main theorem. In Theorem 1, we assume 
that the radius of M is ^ k e for some constant kq > 0. 

Theorem 1. Let Gt denote the time-T map associated with (2.2). Suppose that the function $ defined 
by (2.6) is a Morse function. Then there exist a small constant k\ > and a large constant K2 > n\ such 
that the following holds. If 

(1) I Ai| < m, 

( 2 ) \t\ < K ^> 

(3) ^ > «a, 

then there exists T\ > and a set Ac [T±, 00) of positive Lebesgue measure such that for Te A, Gt admits 
a strange attractor £1 in M and satisfies (SAl)-(SA^) from Section 4- For every interval I C [Ti,oo) of 
length 1, £(A n I) > 0, where £ denotes the Lebesgue measure on K. 

Remark 2.1. The assumption that $ is a Morse function is quite mild and should hold for a typical forcing 
vector field F. We do not formulate precise results of this type in this paper, but such results should hold in 
terms of both topological genericity and prevalence. Prevalence is a measure-theoretic notion of genericity 
that generalizes the concept of 'Lebesgue almost every' to infinite-dimensional spaces. It provides a powerful 
framework for describing generic phenomena in a probabilistic way (see e.g. [6, 7, 10]). 

Remark 2.2. Theorem 1 concludes that Gt exhibits sustained, observable chaos for a set of values of T 
of positive Lebesgue measure rather than for all T € [Ti,oo). This is not a consequence of the nature of 
the proof. Rather, it is a fundamental consequence of the fact that an alternate scenario competes with 
the SRB scenario in the space of T- values. For an open set S of T- values in [T 1; 00), the basin U contains 
a Gy-invariant Cantor set on which Gt is uniformly hyperbolic (a horseshoe) and a periodic sink. The 
trajectory of Lebesgue almost every x £ U converges to the periodic sink. Thus for T £ §, Gt exhibits 
transient chaos: a typical trajectory in the basin will move erratically for some time due to the presence of 
the horseshoe before finally converging to the periodic sink. 



Remark 2.3. The function $ does not depend on the parameters Ai, a, and e. 



Theorem 1 is related to 2 results obtained by Wang and Young in [17]. Wang and Young consider limit 
cycles forced by periodic <5-function kicks. First, they prove that any limit cycle, when suitably kicked, can 
be transformed into a strange attractor. This result is universal but not constructive. An artificially-strong 
kick is needed if geometric conditions are unfavorable for the creation of nonuniform hyperbolicity. Second, 
they prove that the Hopf limit cycle that emerges from a supercritical Hopf bifurcation can be transformed 
into a strange attractor. Here the so-called twist factor plays the role of the shear integral. Unlike the shear 
integral, the twist factor is local in the sense that it depends only on derivatives of the vector field at the 
bifurcation parameter. 

Many of the quantities in Theorem 1 are required to be sufficiently large or sufficiently small. This is an 
unavoidable consequence of the perturbative nature of the analytic techniques used in the proof. However, 
numerical evidence suggests that shear-induced chaos emerges over parameter ranges that far exceed those 
to which the rigorous analysis applies. For example, Lin and Young [9] conduct numerical studies of a linear 
shear flow model previously studied by Zaslavsky [21]. The work of Lin and Young also provides numerical 
evidence that the temporal form of the kicks need not be periodic: temporally-sustained chaotic behavior is 
observed for random kicks at Poisson-distributed times and for continuous-time forcing by white noise. 



3. Derivation of the singular limit 



3.1. Derivation of the normal forms. We derive the normal forms (2.4a)-(2.4b) and (2.5a)-(2.5b) that 
are valid in a small neighborhood of T. For s£§, let {ei(s)}f =1 be an orthonormal basis for E" such that 
e„(s) = j'(s) (where 7' denotes the derivative of 7 with respect to s) and e, is a C 5 function of s for 
all 1 ^ i ^ n. One may choose the first n — 1 vectors in many ways. For example, if 7 is at least C n+5 
and the first n derivatives of 7 are linearly independent, then one may construct the basis by applying the 
Gram-Schmidt procedure to the first n derivatives of 7. For any x G M. n sufficiently close to T, there exist 
unique s € S and y = {yi, . . . , y n _i) such that 



(3.1) 

We use (s, y) as new phase variables. 
Define 



7(s) 



n-l 

E 

»=i 



Viei(s). 



E(s) 



(e 2 (s)) T 



VK(s)) T / 

Differentiating E(s) with respect to s, we have E'(s) = K(s)E(s) where K(s) = (fcj,i(s)) is a skew-symmetric 
matrix of generalized curvatures defined by kj^s) — (e^(s), ei(s)}. If the first n derivatives of 7 are used 
to create E, then this differential equation is the classical Frenet-Serret equation from differential geometry. 
For 1 ^ i ^ n, define the vector 

k 2 ,i( s ) 



ki(s) 



\fc„_ M (s)/ 



Differentiating (3.1) with respect to t, we obtain 
(3.2) 



' JX El? e ^ s ) + § ( VW + E%^.( S ) I =f{x)+eP p>T {t)F{ X ) 



dt 



'-' V j=i 

Taking the inner product of (3.2) with respect to ej(s) for 1 ^ i ^ n — 1 yields 

dy., 



dt 



(f(x),e l ( S ))+eP p , T (t)(F(x),e l ( S )) 



ds 
~dl 



y, hi(s)). 



Taking the inner product of (3.2) with respect to e„(s) yields 
ds 

— ({y, k n (s)) + 1) = (f(x), e n (s)) + eP P;T (t)(F(x), e n (s)). 

Notice that (y, k n (s)) + 1 ^ if is sufficiently small. Consequently, the system 
ds 1 



(3.3a) 
(3.3b) 



((/(sc), e n (s)) + eP PiT (t){F(x), e n (s))) 



dt l + {y,k n {s)) 
^ = (f(x),e l (s))+eP P Mt)(F(x),e t ( s )) - g(y, 



is valid in a small neighborhood of T. 

We now extract the terms of leading order in (3.3a) and (3.3b). For 1 ^ j ^ n, define ipj(s,y) — 
(f(x), ej(s)). For 1 ^ i ^ n — 1, we have 



^( s ,2/) = (^ 1) (.s), 2/ ) + a, y (||y|| 2 ) 



where 



(1) _ d(/(x),e,(s)) 



Here Cs,»(|| J/|| 2 ) denotes a function of s and y for which there exists a constant if > independent of s and 
y such that |0 S ,?/(||2/|| 2 )| ^ ^||2/|| 2 - Expanding tp n {s,y), we have 



where 



Ms, y) = ^ 0) (*) + <^ 1} (*), v) + a, y (||y|| 2 ) 

V4 0) (s) = ll/(7(s))ll 



(1) _ 0(/(aQ,e„(a)) 



H=0 



Set 4>j(s,y) = (F(x), ej(s)) for 1 ^ j ^ n. Writing (3.3a) and (3.3b) in terms of ipj and (pj, when the 
forcing is active (P p ,T{t) = 1) we obtain 



(3.4) 



dt 



ds 



V'n (s) +e4> n (s,y) 
£<f>i{s,y) 



k n (s) 



ds ^\s)+s4> n (s,y) V4 0) (s)+#n(s,y) 



7/>i 0) (s) +e0„(s,y) 
ki{s) ■ y 



y + O s , y (\\yf) 



s<t>i{s, y) 



\ipn , (s) +e<p n {s,y) ) 
When the forcing is off (P Pj t(£) = 0), we have 



(3.5) 



Define 







f/.S 


^° } ( S ) V 








U 0) w 



fc„(s) 



fc„(s) 



^00 



V4 0) (s) +e0„(s,y) 
^00 



2/ + a, y (||y|| 2 ) 



y + o s j\\ y \\ 2 ) 



fci( s ) -y + a^dlyll 2 ) 



6o(«) := 
6i(s) := 



1 

1 



k„{s) 



and let A(s) denote the (n — 1) X (n — 1) matrix with i th row given by 

(3M T - 



In terms of bo, b\, and A, system (3.5) becomes 
(3.6) 



f s =b (s) +<6 1 (a), y) + 0.,y(||»|| 2 ) 



^=A( S )y + a, 9 (||y|| 2 ) 
- as 

Applying the Floquet theorem, there exists a real-valued, periodic (n — 1) x (n — 1) matrix P(s) of period 
2T such that setting z = P _1 (s)j/, we transform (3.6) into 

(3.7a) = 6 °( s ) + (( & i( s )) Tp ( s )) z + M*. z) 

as 

(3.7b) ^=Az + /n(s,z) 

as 

This is the normal form of (2.2) on which we will base our analysis of the flow during the relaxation period 
(when P P: T{t) = 0). We obtain the normal form of (2.2) during the forcing period (when P Pj r(i) = 1) by 
writing (3.4) in (s, z)-coordinates, giving 

(3.8a) £ = b {s) + ((&!(s)) T P(s))z + O s , z (e) + O s , z (ez) + O s , z {\\z\\ 2 ) 

(3.8b) ^ = Az + £P " 1( y ,0) + O s . z {ez) + O s , z {e 2 ) + O s . z {\\zf) 

where 0(s, 0) = (0i(s, 0), . . . , <fi n -i(s, 0)) T . 

3.2. A general form of the singular limit. Let M « T x D be a tubular neighborhood of F in E™, where 
D is a disk of sufficiently small radius so that the normal form (3.8a)-(3.8b) is valid. Let M « T X 4.D. We 
define flow-induced maps iifc : M — ► M and iJ r : M — > M as follows. Let be the time-/? map associated 
with the forced system (3.8a)-(3.8b). We call the 'kick'. Notice that for e sufficiently small, Hf. maps 
M into M. Let H r be the time-(T — p) map associated with the relaxation system (3.7a)-(3.7b). We call 
H r the relaxation map. There exists To = Tq(e) such that if T ^ To, then H r maps M into int(M). The 
composition Gt '■— H r o Hk is the time-T map generated by the flow. Our goal is to show that the family 
{Gt ■ M —> int(Af), T > T } of diffeomorphisms on M has a well-defined singular limit in a certain sense 
as T — > oo. 

Let (so, y ) € M. We write Hk(so, y Q ) = (s, z) and compute H r (s, z). Integrating (3.7b), we have 

z(s) = e {s - s)A (z + J* e-^-^h^r, z{t)) d^J . 

Integrating (3.7a), we have 

f s(T) f s(T) 2 

(3.9) T-p = b (T)dT + z- &i(T) T P(T)e< T -'> A dT + 5^ fc (*(T)), 



fc=l 



where the error terms are given by 

s(T) 

^(s(T)) 



/ biCrfP^e^" 1 )* / e" (£,*(£))d£ or 

£ 2 (s(T)) = / /lafr.zMJdT. 



Letting T — )• co in (3.9) yields nothing meaningful. However, we use the fact that s is can be computed 
modulo 2T to introduce an auxiliary parameter a € S and thereby obtain the singular limit. Recall that po 
is the period of r). As a varies from to 2T, 7 traverses L 2 times. Let t : [0, 2L) — > [0, 2po) be the strictly 

7 



increasing function denned by r)(t(a)) = 7(a). For m E Z + and a <G S, set T — p + 2pom + t(a). Substituting 
into (3.9), writing s{p + 2pom + t(a)) = s + 2Lm + s(p + 2pom + i(a)), and using the fact that 

D+2Lm 

60 (t) dr = 2pQin 

J V 

for all v E R, we obtain 

ps+s(p+2p m+i{a)) /• s(p+2p m+i(a)) 

i(a)= b a (T)dT + z- fe 1 (r) T P(r)e (r ~" )A dT 



(3.10) 



^2E k (s(p + 2p m + i(a))- 



k=l 



Define G a ^ m -i : M — > int(M) by G a m -i(s , y ) = (s(p+2p m+t(a)), y(p+2p m+i(a)). It follows from [17, 
Proposition 3.1] that there exists Soo(soj ?/o> a ) such that 

lim s + s(p + 2p m + i(a)) = s^So, y , a) 

m—>oo 

and Soo(so, Vo,a) is defined implicitly by taking the m — > 00 limit in (3.10): 

/■Soo(SO,!/0: a ) If 00 \ J^-, 

(3.11) *(«) = / b (T)dT+^z,j^ b 1 (r) T P(T)e<r->»dT}+Y l Ek(oo). 

The family of maps {G a! o : M — >• T x {0}} ae s defined by 

G a ,o(s , Vo) = («oo(so, 2/ , a )> °) 

is the desired singular limit. It follows from [17, Proposition 3.1] that the maps 

(so, 2/ , a ) !-> G a , m -i (s , |/o) 

converge to the map 

(s , 2/o, a ) | -> G a ,o(s ,2/o) 

in C 3 (M x S) as m -> 00. 

3.3. A computable form of the singular limit. ^From this point forward, we assume the setting of 
Theorem 1. We now extract the primary terms in the right side of (3.11). Recall that the shear integral £ 
is defined by 

I-2L 

\ dr. 



S = (S 1 ,...,S„_ 1 )= / 6i(r) T P(r). 
Jo 



and that the shear factor is given by a = ||S||. We assume that the operator A is diagonalizable and that 
the z-coordinate has been chosen such that A = diag(Ai, . . . , A n _i), where > Ai ^ A2 ^ ■ • • ^ A„_i are 
the eigenvalues of A. Fix the normalized shear vector — and the eigenvalue ratios p,i — 4^ for 1 ^ i ^ n — 1. 
Set p = 1 for notational simplicity. We regard tr, e, and Ai as the parameters associated with the singular 
limit. 

Expanding the second term on the right side of (3.11), we have 

/>oo />oo />oo 

/ 6 1 (r) T P(r)e( r -^ A dr = / Xe^-^ dr + (6!(t) t P(t) - S)e (r ~ s)A dr 

J S J S J § 

= ~d + / 00 (b 1 (T) T P(r) - £)e< r -^ A dr, 

J s 



(3.12) 



where 

d= r^-^dr= (-^ 



» / i=l 



Let Hk '■ M —> M be the time-1 map generated by the system 
(3.13a) ^ - b (s) 



(3.13b) 



ds 

dz _ eP^ 1 (s)^(s, 0) 
ds ~ Vi 0) (s) 



obtained from (3.8a)-(3.8b) by retaining only the terms of leading order. For (so, Vq) & M, write Hk(so, y ) 
(s,z). Integrating (3.13a) and (3.13b) gives 



(3.14) 



1= / 6 (r)dr, 

J So 

r P- 1 (r)0(r,O) 

z — zq -f- e 



dr. 



^>0 (r) 

Proposition 3.1. There exists a system constant Kq > such that 
(3.15) s = s + ^i(s ,y ), z = z + £ 2 ( s o, Vo) 

w/iere 

ll&KVo = 0}llc3(s) < ^o£, ||€ 2 |{tfo = 0}llc3 ( s) < ^oe|Ai|. 
Setting y = 0, define g(so,a) — Soo(so,0,a). Substituting (3.12), (3.14), and (3.15) into (3.11), the 
value g(so, a) is defined implicitly by 



t(a) + 1 = 



(3.16) 



6 (r) dT+((z+| 2 (s o ,0)),d) 



/ b (r)dT + z- (b 1 (r) T P(r)-S) e ( T -^ A dr + ^S fc (oo) 
^ fe=l 



Rescaling d, we define 



giving 



n-l 



so V% (t) 



dr ) , 



The higher-order terms are given by £i = £1(00), £2 = £2(00), 

£ 3 = ^°°(Mt) t P(t) - £)e (r ^ )A dr^ , £ 4 = - f 6o(t) dr, £ 5 = (£ 2 (s , 0), 2). 
Setting £ = ^fc=i &k and substituting into (3.16), we obtain the final form of the singular limit: 

bo(T)dT+-—<P( So ) + Z. 



(3.17) 



t(a) + 1 



|Ai, 

Proposition 3.2. There exists a system constant K\ > such that the following hold. 



ae ( e 
ere ( e\ 

£ 2 ||C3(S)<^1^(-J 



as 



£3||c»(s)<^in-r(|Ai|) 

|Ai| 



£4||c 3 (S) ^ Ki 



ae (\\ 



|Ar| V 



as 



£ 5 ||c3 ( s)<^ip^(|Ai|) 



4. Theory of rank one attractors 



Let D denote the closed unit disk in 1" and let M = S x D. We consider a family of maps G a 



i, 



M M, where a = (a\, . . . , a^) E V is a vector of parameters and b E Bq is a scalar parameter. Here 
V = Vi x • • • x Vj. C R fc is a product of intervals and B C K \ {0} is a subset of R with an accumulation 
point at 0. Points in M are denoted by (x, y) with x E S 1 and y E D. Rank one theory postulates the 
following. 

(HI) Regularity conditions. 

(a) For each b g £?o, the function (x,y, a) n- G a .b(x,y) is C 3 . 

(b) Each map G a .b is an embedding of M into itself. 

(c) There exists > independent of a and 6 such that for all a E V, b E B , and z, z' E M, we 
have 

ldet£>G a , b (z)| 
|detDG„ lt (a/)| " D ' 

(H2) Existence of a singular limit. For a 6 V, there exists a map G a ,o '■ M — > S 1 x {0} such that the 
following holds. For every (x, y) E M and a E V, we have 

lim G a b (x, y) = G a (x, y) 

6 — s-0 

Identifying S 1 x {0} with S\ we refer to G a fl and the restriction f a : S 1 — >• S 1 dehned by /o(x) = 
G„ o(^|0) as the singular limit of G a j,. 

(H3) C 3 convergence to the singular limit. We select a special index j E {1, . . . , k}. Fix a, € Vj for 
i 7^ j. For every such choice of parameters dj, the maps (x,y,aj) >-> G Ut i,{x,y) converge in the C 3 
topology to (x, y, a,j) H> G a ,o{x, y) on M x Vj as & — > 0. 

(H4) Existence of a sufficiently expanding map within the singular limit. There exists a* = 
(a*, . . . , at) E V such that /„* E M, where M is the set of Misiurewicz-type maps defined in Defini- 
tion 4.1 below. 

(H5) Parameter transversality. Let C a * denote the critical set of / a ». For a 3 E Vj, define the vector 
a,j E V by hj = (a*, . . . , a*_ 1 ,aj, a* +1 , . . . , a* k ). We say that the family {f a } satisfies the parameter 
transversality condition with respect to parameter a, if the following holds. For each x E C a * , 
let p — f a *(x) and let x(a,j) and p{ctj) denote the continuations of x and p, respectively, as the 
parameter aj varies around a*. The point p(a,j) is the unique point such that p(a,j) and p have 
identical symbolic itineraries under /„ and / a * , respectively. We have 

d_ 

daj 



daj 



3 



(H6) Nondegeneracy at 'turns'. For each x E C a *, there exists 1 ^ m ^ n — 1 such that 

d 



-G a . fi {x,y) 

oy m 



7^0. 

y=0 



(In 



(H7) Conditions for mixing. 

i, 

(a) We have e3 > 2, where Ao is defined within Definition 4.1. 

(b) Let Ji, . . . ,J r be the intervals of monotonicity of /„« . Let Q — (qi, n ) be the matrix of 'allowed 
transitions' defined by 

fl, if f a *(Ji) => J m , 

1 0, otherwise. 
There exists N > such that > 0. 
We now define the family M. 

Definition 4.1. We say that / E C 2 (S 1 ,]R) is a Misiurewicz map and we write / E M if the following hold 
for some neighborhood U of the critical set C = C(f) = {x E S 1 : f'(x) = 0}. 
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(A) (Outside of U) There exist A > 0, M € Z+, and < d sC 1 such that 

(1) for all m ^ M 0) if .f (x) £ U for < i < m - 1, then |(/ m )'(x)| ^ e A ° m , 

(2) for any m € Z+, if (x) ^ C/ for sC m - 1 and f rn (x) e 17, then |(/ m )'(x)| > d e A ° m . 

(B) (Critical orbits) For all c G C and i > 0, p(c) <£ U. 

(C) (Inside U) 

(1) We have /"(x) ^ for all x G U, and 

(2) for all x € f7 \ C, there exists p (» > such that f (x) g [/ for all i < p («) an d |(/ Po(:r) )'(x)| > 



d -l e 3 A oPo(x)_ 

Rank one theory states that given a family {G 0i f,} satisfying (HI)— (H6), a measure-theoretically signifi- 
cant subset of this family consists of maps admitting attractors with strong chaotic and stochastic properties. 
We formulate the precise results and we then describe the properties that the attractors possess. 

Theorem 4.2 ([15, 18]). Suppose the family {G a .b} satisfies (HI), (H2), (H4), and (H6). The following 
holds for all 1 j ^ k such that the parameter aj satisfies (H3) and (H5). For all sufficiently small b € Bq, 
there exists a subset Aj C Vj of positive Lebesgue measure such that for aj G Aj, Ga^t admits a strange 
attractor n with properties (SA1), (SA2), and (SA3). 

Theorem 4.3 ([15, 16, 18]). In the sense of Theorem 1^.2, 

(H1)-(H7) (SA1)-(SA4). 

Remark 4.4. The proof of Theorem 4.2 for the special case n = 2 appears in [15]. The additional com- 
ponent (H7) =>■ (SA4) in Theorem 4.3 is proved in [16]. For general n, Wang and Young [18] prove the 
existence of an SRB measure for Ga<,6 if a j € Aj. The complete proofs of (SA1)-(SA3) (and (SA4) 
assuming (H7)) for Ga,^;, with aj € Aj will appear in [14] for general n. 

We now describe (SA1)-(SA4) precisely. Write G — ,b- 

(SA1) Positive Lyapunov exponent. Let U denote the basin of attraction of the attractor f2. This 
means that U is an open set satisfying G(U) C U and 



O = p| T m (U). 



For almost every z £ U with respect to Lebesgue measure, the orbit of z has a positive Lyapunov 
exponent. That is, 

lim — log||iX? m (z)|| > 0. 
ro-»oo m 

(SA2) Existence of SRB measures and basin property. 

(a) The map G admits at least one and at most finitely many ergodic SRB measures each one of 
which has no zero Lyapunov exponents. Let 1/%, • • • ,v T denote these measures. 

(b) For Lebesgue-a.e. z € U, there exists j(z) £ {1, . . . , r} such that for every continuous function 

tp : U -> R, 

- m— 1 „ 

m * — ' / 

(SA3) Statistical properties of dynamical observations. 

(a) For every ergodic SRB measure v and every Holder continuous function tp. : — > R, the 
sequence {tp o G l : i G Z + } obeys a central limit theorem. That is, if J tp dv — 0, then the 
sequence 

m— 1 



i=0 

converges in distribution (with respect to v) to the normal distribution. The variance of the 
limiting normal distribution is strictly positive unless tp = ip o G — tp for some ip £ L?(y). 

li 



(b) Suppose that for some N ^ 1, G N has an SRB measure v that is mixing. Then given a Holder 
exponent 77, there exists t = r(^) < 1 such that for all Holder ip, i/j : fi — > K with Holder 
exponent 77, there exists if = if (<y9, ?/;) such that for all m £ N, 



(ipaG mN )il>dv- I (pdu I ^dv 



(SA4) Uniqueness of SRB measures and ergodic properties. 

(a) The map G admits a unique (and therefore ergodic) SRB measure v, and 

(b) the dynamical system (G, v) is mixing, or, equivalently, isomorphic to a Bernoulli shift. 

5. Verification of the rank one hypotheses 

We view the singular limit {G a .o : a £ S} as a function of 3 parameters: e, er, and Ai. We show that the 
family {G a m -i : a £ S, m £ Z + } satisfies (H1)-(H7) if the parameters e, <r, and Ai satisfy certain scaling 
assumptions. 

5.1. ID analysis: verification of (H4), (H5), and (H7). Recall that g(s,a) is defined implicitly by 



t(a) + l = J &o(r)dr+|— *(«) + £. 

Defining / Q (s) = g(s, a), A = 4^-, and \&(s) = $(s) + A _1 £, the singular limit becomes 

rfa(s) 

(5.1) t(o) + l= / 6 (r)dr + Atf(s). 



For a map / : § -> § and 6 > 0, let C(/) = {s : /' (s) = 0} and let C s (f) = {s : \s-s\ < 6 for some s £ C(/)}. 
We assume the following about ^: there exist positive constants if 2, do, d±, and g?2, and a constant <5o 
satisfying < do < \d\, such that the following hold. 

(Al) ||*||c«(s) < #2 

(A2) >do for sGC 5o (*) 

(A3) If *'( Sl ) = #'(s 2 ) = and s 1 f s 2 , then \s x - s 2 \ > d x . 
(A4) \V{s)\ >d 2 far«€S\C tfo (#) 

Because $ is a Morse function, Proposition 3.2 implies that assumptions (A1)-(A4) are satisfied if a ^ 1, 
I Ai| is sufficiently small, and r^-r is sufficiently small. 

We now compare the map f a to the map 4". Let {v\, . . . ,v qo } be the set of critical points of 'J. Set 

6 = 

Lemma 5.1. There exists Aq > and positive constants if 3, K4,, and if 5 such that the following hold for 
fixed A > Ao . 

(a) C{f a ) = {vi, . . -,v qo } with \v t - v { \ < if^A -1 for 1 < i q 

(b) \f!f(s)\>K i A for alls G Q(/ a ) 

fcj |/»| > K 5 Ai for all s £ 3 \ G |f (/ a ) 

Proof of Lemma 5.1. Differentiating (5.1) with respect to s, we obtain 

(5.2) b (s) - A*'(s) = b (f a (s))f' a (s). 

Setting f' a (s) = gives bo(s) = A^'(s). Since 60 is bounded above and bounded away from 0, (A2)-(A4) 
imply (a). Solving for /„(s), we have 

(5.3) m - M f'r f f» w 

On §\ Ci { (/„) we have |*'(a)| > if£ using (a), (A2), and (A4). Estimate (c) now follows from (5.3). 
Differentiating (5.2) with respect to s, we obtain 

(5-4) b' ( S ) - A*» - &' (/ a (s))[/^)] 2 = boUaisMis). 



For all s £ Q(/ ), we have |*'(s)| < A£ by (Al) and (a). This implies that |/£(s)| < AAi on Q using (5.3). 
Therefore the second term on the left side of (5.4) dominates and (b) holds. ■ 

5.1.1. Critical curves. Assume A > A and let A C § be a parameter interval. For a £ A, we have 
C(/a) = {vi(a), . . . ,v qo (a)} by Lemma 5.1. Write 7^(a) = Wj(a) for 1 ^ i < g . For 1 ^ fc ^ g and ieN, 
define 7 f ) (a):=/:(7W(a)). 

Differentiating j[ (a) = f a (j^ (a)) = f(^ k H a )i a ) with respect to a, we have 



(5.5) 



^W = g(^W,a).A 7 ( fc ) (a) + g (7W(a) , a) 



a/ 

9a 



(7 W (a), a). 



Differentiating (5.1) with respect to a and using the fact that £i(a) — 60(a), we obtain 
(5.6) 
Thus 



— f( \ = b °( Q ) 
da J{S ' 1 b (f( S) a)Y 



d ( k )t \ ^ min s6S &o(s) ^ n 

Til w > rr\ > °- 

da maxjgg o (s) 

More generally, an estimate on j-7|S(a) for i £N follows from the recursive formula 



(5.7) 



g( 7 <" w ,«,.^(„, + |( 7 f(«), o) . 



Lemma 5.2 (Growth estimate for derivatives of critical curves). There exists Ai ^ Ao such that the following 
holds for all A > Ai . For any k 6 {1, . . . , <7o} an d * £ N such that 7 j (a) € S \ Cf (\&) /or all 1 ^ j ^ i, then 

d 



(5.8) 



( fc ) / 1 



da 



K 5 k 1 

I f A, 



> As . 



Proof of Lemma 5.2. Estimate (5.8) follows from (5.7), estimate (c) from Lemma 5.1, and the fact that for 
all sgS and a £ A we have 

d , _ max seS 6 (s) 



7(s,a) < 
da mm seS o (s) 



A fi 



Lemma 5.3 (Distortion estimate for critical curves). There exists A 2 Ai and > swc/i i/iat tte 
following holds for all A > A 2 . For any k £ {1, . . . , go} a77 ^ ^ 2, fei A 6e a parameter interval such 
that 

(a) j\ k) (A) C S \ Q(*) /or 1 < i < ra - 1, and 

f&j ^(7^ 1 (A)) < £ (£ denotes Lebesgue measure on E>). 

Then for all a, a £ A, we have 



(5.9) 



da '« W 



da' n w 



<D t . 



If n = 1, £/ien (5.9) ZioWs /or all k £ {1, . . . , go} and for all a, a £ §. 
Proof of Lemma 5.3. For n = 1 and a, a £ S, the estimate 



ibi fe) («) 



< A fi 2 



follows from (5.5) and (5.6). For n ^ 2 and a, a G A, let = 7^* (a) and s» = 7^(0.). We have 







da** 








da* 1 





fL(Si-l) ■ £si-l + £/a(Si-l) 



1 + e>(A 
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This implies the estimate 



log 


^s 


= log 


da 1 




da bn 




da 1 






^ 1 

1=1 1 



Y, log 



f'M 



+ 0(1) 



5^1og(l + 0(A" 



The equality 

\f'a(Si)-f'am = 



1 



((6o(si+i) - 6o(si+i))/a(si) + A(*'(«i) - + (6 (si) - 60(50)) 



implies the estimate 



log 



cia r 



i—1 i—1 i—1 

/ ri— 1 \ n— 1 n— 1 

AT |s„ - s n | + \si - 8i\ + KAi l*< " «i| + ^ A ~* 12 \ § * ~ s *\+ 



ra-3 



n-2 



< A>„ - S„| + A> n _! - s n _i| 2 A-i* + K\a n -! - S„_i|(A5 + A"*) £ A"i' + 0(1) 

8=0 1=0 

= 0(1). 



5.1.2. Verification of (H4): Definition 4.1(B). We prove the existence of a parameter a* such that f a * 
satisfies Definition 4.1(B). We will then show that if A is sufficiently large, then for any parameter a, if f a 
satisfies Definition 4.1(B), then f a 6 M. 

Proposition 5.4. There exists A3 A2 such that if A ^ A3 and A C S is a parameter interval satisfying 
t(A) = 3D 1 K 6 q £, then there exists a* £ A such that for all c G C(/ Q .), (c) G § \ Q(*) /or aZZ neN. 

Proof of Proposition 5.4- We inductively construct a nested sequence of parameter intervals A = Ao D Ai D 
A 2 D • • • such that a* G f)i=o ^« nas * ne desired property. 



Definition 5.5. The (go + l)-tuple (A n ;ii yJl ,. 



i) is called an admissible configuration if A„ is a 



subinterval of Ao and if for every k G {1, . . . , go}, *fc,n ^ n and the following conditions are satisfied. 

(Ml) 7 f } (A„) n C c (*) = for all t < * fcjn 

(M2) For all a, a G A n , we have the distortion estimate 

^--/ (fe) fa) 



^T (fc) fa) 



<Di. 



(M3) £( 7 W +1 (A n ))>3A^ 

We inductively construct admissible configurations for all n G N such that ik lU — > °o as n — > 00 for every 
k. We begin with n = 1. Let 

d := min Is — £|. 

s,teC(*) 

We assume that 3-Di-ftT|g £ < 2^- -^ et *M = ^ or au ^- We d 100136 Ai as follows. We have 

(fe), s _ 60(a) 

^"^(^(a))' 



Consequently, 7} (Ao) meets at most one component of C^(^) and we have 

^((^lAo)- 1 ^^))) 2 
*(A ) "3g " 

Even in the worst-case scenario in which the g intervals {(71 ) _1 (Cf(^ r )) : 1 ^ fc ^ go} are evenly spaced 
in A , there exists a subinterval Ai of A with £(Ai) ^ Dl g f+i" g such that 7^ ) (A i ) n = for all fc. 

Property (Ml) holds by design and (M2) follows from Lemma 5.3. Property (M3) holds if A is such that 

£( 7 f ) (A 1 ))^A^(A 1 )^3^ 1 g e 

Now assume that for n € N we are given an admissible configuration (A„; ii n , . . . , iq n)- We construct 
an admissible configuration at step n + 1 as follows. Partition the set {1, . . . , go} into 2 sets: A, the set of 
indices that are 'ready to advance', and {1, . . . , qo}\ A, the set of indices that are not ready to advance. The 
index fc is in A if (II) and (12) hold: 

(11) (■(ritk (An)) < £ (distortion estimate holds for the next iterate) 

(12) li^i^ j_i(A n )) < \d (image of the next iterate meets at most one component of C^(*)) 

Suppose that A ^ 0. In this case, set 

. = Uk,n + l, ifke A; 
' n+ [ik.n, if k E {1, . . . , g } \ A. 

We now find A n+1 so that (A„ +1 ; ii,„+i, . . . , i qo . n +i) is an admissible configuration. Let fc € A. Using (M3) 
and (12), we have 

3D iqo ^e(^l +1 (A n ))<^d. 

This implies that the fraction of 7^ +1 (A„) in C^^) is bounded above by 3 ^ qo ^ = 3D iqg ■ Using (II) and 
Lemma 5.3, we have 

£(( 7 g ) „ +1 |A»)- 1 (Q(^))) 2 
£(A n ) ^ 3g ' 

Arguing as in the n — 1 case, there exists a subinterval A„ +1 of A n such that £(A n+ i) 3 ^ +1 ^ £(A n ) and 

for all fc e A, ^\ +1 (A n+1 ) n Q(^) = 0. For k e A, (Ml) holds by design and (M2) follows from (II). The 
inequality 

implies that (M3) holds if A is such that 

%£ +2 (A„ +1 )) > (^) (f A*) > 3D iqo Z. 

Now let fc G {1, . . . , go} \ A. Properties (Ml) and (M2) are inherited from step n. If (II) fails for index fc, 
then (M2) gives 

so index fc satisfies (M3) if A is such that 

«i"„ +1 (A„ + ,)) > (jjs^j) (f Ai) > 3D, 9 „f. 
If (II) holds but (12) fails, then (M3) holds for index fc if A is such that 

<1 +1 (A„ +1 )) > (^pTj) (^) > 3A, e 
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If A = 0, then let A' n be the left half of A„. We claim that (A^; i x . iq a ,n) is an admissible config- 

uration. For each index k, properties (Ml) and (M2) trivially hold. Property (M3) is established (for A 
sufficiently large) by arguing as above in the 2 cases 

(1) (II) does not hold, and 

(2) (II) holds but (12) fails. 

Repeat the halving process until A ^ 0. ■ 

5.1.3. Verification of (H4): f a * satisfies Definition 4.1(B) f a * G M. We show that if A is sufficiently 
large and a* £ § is as in Proposition 5.4, then f a * £ M. This implication is a consequence of Lemma 5.1 
and the following binding estimate. 

Proposition 5.6. There exists Ki > such that for A sufficiently large and a* as in Proposition 5.4, we 
have the following. For c £ C(/ a ») and s £ § satisfying \s — c| ^ A~T2 ; l e t rn{s) be the smallest value of 
m£Z+ such that \f£(s) - /™(c)| > |£. Then m{s) > 1 and 

\(C {S) )'( S )\^(K 7 A)^. 
Proof of Proposition 5. 6. We begin with a spatial distortion lemma. 

Lemma 5.7 (Spatial distortion estimate). There exists D2 1 such that the following holds for all a £ S. 
For s,s £ S, let m £ Z + be such that 7Tj, the segment between fl(s) and f l a {s), satisfies £(iti) < j£ and 
■Ki H CiAfa) — for all ^ i < m. Then 



(/a" l )'(s) 

Proof of Lemma 5.7. Writing = /„(s) and Sj = /„(s) and using Lemma 5.1 and its proof, we have 



log 



(/« m )'00 



m—1 

E 

i=0 



log 



E 

i=0 



< i^A"' ( AV 



A|l^ 7 Hc«(s) 
minxes b (to) 



m—1 

E 

i=0 



< (ifA~ 
= 0(1). 



m— 1 

iTA3)|s m _i - s m _i| ^(K 5 AJ)" 

i=0 



Returning to the proof of Proposition 5.6, write / = /„*. We first show that m(s) > 1. We have 

l/(s)-/(c)|^|/"(C)l(s- C ) 2 
for some £ satisfying |£ — c| ^ A~n. Arguing as in the proof of Lemma 5.1, |/"(C)| ^ ATA. Therefore 

|/( S )-/(c)KiTA-f <| 

for A sufficiently large. 

Now assume m(s) > 1. Using Lemma 5.7, we have 



< |/ mw (s)-/ m ^(c)| 

= l(/ m(sK1 )'(Ci)l • |/(s) - /(c)| (for some Ci between /(s) and /(c)) 
^^ 2 | (/ ^)-iy (/(c)) |.| /(s) _ /(c) | 
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and therefore 

(5.10) e < D 2 \(r^y(f(c))\ ■ i /"(o i ■ (s - c) 2 . 

Reversing inequality (5.10) at time m(s) — 1, we have 

(5-11) £ > D^\{r {s) - 2 )'(J{c))\ • | /"(C) I ■ (* - c) 

Estimating |(/ m(s)_1 )'(/(c))| from below using (5.10) gives 
|(/™W)'( S )| = |/'( S )-/'(c)|-|(/ m W- 1 )'(/( S ))| 

^ £>2 ^(/^"^'(/(c))! • |/"(C 4 )| • I* - c| (for some d between s and c) 



(5-12) ^ n2 

i° — m \ u 

l/"(q)l 



D 2 2 \s-c\ VI/"(C)| 



Arguing as in the proof of Lemma 5.1, Tji/Wj\ ^ A' > since C4 and C are between s and c. Using this fact 
and estimating |s — cj^ 1 from below using (5.11), (5.12) implies 



Kr wy (s) | ^ ^ / ^-^(/"^-^-(/(c))!-!/"^)! 



^ (ata^-s 

> (ATA)t^. 



5.1.4. Verification of (H5) and (H7). The following lemma facilitates the verification of (H5). 
Lemma 5.8 ([12, 11]). Let f = /„*. Suppose that for all x 6 C(/ a »), we have 



fc=0 

T/ien /or each x £ C(/ a *), 

f> [(9a/a)(/ fc (^))]a^ 
fr (/*)'(/(*)) 



^ l(/ fc )'(/(^))l <0 °' 



^-fa{x{a)) - ^-p{a) 
da da 



Property (H5) follows from Lemma 5.8 for A sufficiently large. To see this, suppose f a * g M and let 
c G C(f a *)- For k e Z + , we have 

K/aO'(/(c))|^(^ 5 Ai) fc 

by Lemma 5.1(c). Since K^ 1 ^ ^p/(s, a ) ^ -^6, we conclude that if A is sufficiently large, then 

fro tenfa-w - K ° -h(K,Ai) k 

Property (H7) follows from Lemma 5.1 and Proposition 5.6 provided A is sufficiently large. 

Appendix A. Some proofs 

We assume throughout Section A that L = 1. Notice that if V denotes a vector field, then 

= ^ 1 fl!||z|| 2 = z = _z_ 

1 ' j tfc tfo 2||z|| ||z|| ' ds \\z\\ ' 

We will use this fact together with the following Gronwall-type inequality: 
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Lemma A.l. Assume that j3 is a constant, the function (p is continuous on the interval [s,s], and that the 
function u is differ entiable and satisfies ^ ^ j3u + <p on (s. s). Then, for all s G (s, s), 

u{s) ^ u(s)e 0{s -~ s) + f e 0{s - T) ip{T)d T . 



Proof. Suppose v(s) = it(i)e^ s s ~> + f-e 13 ^ T ^tp(T)dT. Then v satisfies the equation fj(s) = fiv(s) + ip(s) 
with v(s) = u(s). Since u — v is differcntiable, -^(u — v) ^ f3(u — v), and (u — v) (s) — 0, a standard Gronwall 
argument shows that u ^ v. I 

We get immediately 

Corollary A. 2. Suppose that in Lemma A.l ^g(s) sC + C a e Xl ^- s \ Then 

««<(«(«) + ^)e*<«-'>. 

Our first application of a Gronwall inequality is 

Lemma A. 3. Assume z solves the forced equation (3.8b) with z(sq) — Zq and fix a constant K > 0. // 
e/|Ai| is sufficiently small, 

(A.2) \\d?d l So z(s)\\^Ce (0<Z + m<3) 

as long as s — Sq ^ K . Moreover, 

(A.3) IKz-llKClArl. 

Proof. Equation (3.8b) reads 
dz 

(A.4) = Az + h 3 (s,z) with h 3 (s,z) = O s (e) + O s , z (ez)+O s , z (£ 2 )+O s , z (\\z\\ 2 )- 

ds 

Assuming ||z||/|Ai| and e/|Ai| are sufficiently small, (A.l) implies 

d W z W < A i,| r || , r . 
— <: Y \\z\\+C e. 

By Lemma A.l, 

||z(,s)|| < ||2 ||e^ (s - So) + (l - e^K-'")) sC \\z Q \\eM^-^ + C s(s - s Q ). 

For s — s ^ K we get ||z(s)||/|Ai| ^ ||zoll/l^i| + C Ke/\\i\, which proves the assumption legitimate. 

Differentiating (A.4) with respect to s up to two times yields an expression for d™z(s). One immediately 
obtains 

\\d™z(s)\\^Ce 

for < m < 3 and s — s ^ K. 
Equation (A.4) implies 

(A.5) z(s) = e^-^zo + f e^ s - T ^h 3 {T,z{T))dT. 



Differentiating this with respect to So up to three times and evaluating at s = sq yields 
d so z(s ) = -Az - h 3 (s , z ) 

d 2 Q z(s ) = A 2 z + Ah 3 (s , z ) - d s h 3 (s , z ) - Dh 3 (s , z )d So z(s ) 

d 3 So z{s Q ) = -A 3 z Q - A 2 h 3 {s a , z Q ) + 2Ad s h 3 (s , z ) - d 2 h 3 {s , z ) + ADh 3 (s , z )d So z(s ) 

- D(d s h 3 )(s , z )d So z(s ) - Dh 3 (s , z a )J-d So z(s ) - D 2 h 3 (s , z )(d So z(s Q ), d So z(s )) 

aso 

- Dh 3 (s ,z Q )d 2 z(s ). 
Clearly, for 1 I 3, 

||# z(*o)|| < Ce. 
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Such initial conditions are needed for analyzing the variational equations 

(A.6) ~d S0 z = (A + Dh 3 (s,z))d S0 z 
as 

(A.7) ^d 2 So z = (A + Dh 3 (s, z)) d 2 So z + D 2 h 3 {s, z)(d Sa z,d SQ z) 

^d^z = (A + Dh 3 (s, z)) d^z + 3D 2 h 3 (s, z)(d 2 Q z,d So z) + D 3 h 3 (s, z){d So z,d SQ z, d Sa z). 
One then checks recursively, using (A.l) and Corollary A. 2, that 

||fli o *0»)|| < Cee^ {s - So] ^ Ce 

hold for 1 ^ I $J 3 and s — sq ^ K. 

Equations (A.6) and (A.7) provide us with an expression for d s d l s z(s) with I = 1 and I = 2. Moreover, 
(A.6) can be differentiated with respect to s to yield an expression for d 2 d So z(s). The bounds in (A. 2) are 
then readily obtained. 

Finally, we will prove (A. 3). To this end, notice that 

(A.8) ~d Zo z = (A + Dh 3 (s, z)) d Zo z. 

In particular, each row, d Zoi z, of d ZQ z satisfies this equation. Hence, by principle (A.f), ^||<9 Zoi z|| ^ 

4p||9 Zf) t z\\, so that the matrix d Zo z remains perpetually bounded. Integrating both sides of (A.8) from So 
to s and recalling d Zo z(s ) — 1 gives 

||0, a z(s)-l||<|s-a o |(|A n -i|+ sup \\Dh 3 ( s ',z)\\ sup \\d Zo z(s r )\\) , 

where || • || now denotes the matrix norm induced by the Euclidean norm. This estimate implies (A. 3). I 

Proof of Proposition 3.1. Throughout the proof, || • \\c$ will stand for the C 3 -norm with respect to sq. 
By (3.8a) and (3.14), s and s have to satisfy 

(A.9) f b Q (T)dT = p= f b Q (r) + v(t) dr, 

where v(s) = b]{s)P{s)z(s) + O s . z (e) + O s , z {e\\z{s)\\) + O s . z {\\z{s)\\ 2 ) and z = z(s) solves (3.8b) with 
z(s ) = z . 

We use the implicit function theorem to find s. Clearly, F:KxI->R: (so,s)4 J &o(t) dr — p is 
C 3 . Observe that F(sq, sq) — —p and hm s _i. 00 F(so, s) — oo as minfeo = m > 0. By the intermediate value 
theorem, there exists a number s such that F(sq,s) = 0. Because J^F(so,s) = bo(s) J? m, the implicit 
function theorem implies that s is a C 3 -function of sq. Notice that F(so + l,s + 1) = F(sq,s), so that 
s(sq + 1) = s(sq) + 1 which implies that sq i-» s(sq) — sq is periodic. 

Now that we have s, let us define the function 

g(0 :=-p+ / 6q(t) + v(r) dr. 

Notice that, denoting £i = s — s, the right side of (A.9) is equivalent to <?(£i) = 0. The Taylor expansion 
9(0= 9(0) +9 , (0)S + 6 3 g(0 yields 

G(0-=- g L,(9(0) + S 2 9(0)=^ 

which we regard, for all fixed Zq, as a fixed point equation on the space of C 3 functions £ = £(so). Assuming 
G is a contraction in a closed, origin-centered, ball B r C C 3 of radius r, there exists a unique solution, £i, 
to G(£) = £ inside the ball. Next, we prove that for a suitably small value of r, G is indeed a contraction. 
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First, notice that 



v(t) dr = (s — sq) / v((l — t)sq + ts) dr 



9(0) 



g'(0) = b (s)+v(s) 
$29(0 = e \ (1 - r) s"(£r) dr = £ 2 / (1 - r) (b' + v')(§ + £r) 



are smooth functions of so- Because s is C 3 in s and inf SQ </(0) > 0, the bounds (A. 2) yield 



ff'(0) 



s$ C and 



7(0) 



ff'(0) 



< Ce. 



c 3 



Moreover, 



IIMOIIc < Cll^llca sup ||(6 + */)(§ + C)I|C3. 



Hence, ||G(£)||c 3 ^ Co(e + r 2 ) for some Co- Choosing r = 2Co£, we have G(B r ) C -B r for e small enough. 

Second, let £ and £ 2 be elements of Since the map £ t— > G(£) is differentiable and the operator norm 
of the derivative obeys the bound supf & § r \\DG(0\\z(C 3 ) ^ C su P^es r ll^^2ff(C)IU(c 3 ) ^ Gr, the mean value 
theorem yields HG^ 1 ) — G(£ 2 )|jc3 ^ Gr||£ 1 — £ 2 ||c 3 - Hence, G is a contraction on i? r if e is sufficiently small. 

We will now prove that the fixed point, £i, of G is a periodic function of sq- Let us denote z(s, Sq, Zq) 
the solution and w(s)| so the function v defined above, when the initial condition z(sq) = Zq is being used. 
Because h 3 (s + 2, z) = h%(s, z) in (A. 4), we have z(s + 2, s + 2, z ) = z(s, s , z ) and u(s+2)| So+ 2 = v(s)\ ao . 
Since g(£±) — for all values of s and s(s + 2) = s(sq) + 2, the computation 



fl(&)(*o + 2) = -p + 



s(s +2)+£i(s +2) 

s +2 

(*o)+£i(«o+2) 



& o(t) + w(r)| Sn+2 rfr = -p + 



s(s )+2+4i(s +2) 



^o(t) + u(T)| So+2 dT 



s +2 



•s(&)(*o) + 



6 (t + 2)+i;(t + 2)| So+ 2 rfr 

s(so)+£i(s +2) 



5(so)+Ci(so+2) 



M r ) + w(r)| So dr 



&o( r ) +w(r)| So dr, 

's(so)+fi( s o) 

implies that the last integral vanishes despite the fact that the integrand is positive, so we must have 
£i(so + 2) =6(*o). 

As the last step, we will bound the difference z — z. Let z^ 1 ' and z^> solve (3.8b) and (3.13b), respectively, 
with the initial condition z^(sq) — z^(s a ) = z a . Both of these are C 3 functions of (s ,Zq) by the 
smoothness of the vector fields. By definition, z = jjWfs) and z = z^ 2 \s). We need a bound on the G 3 
norm of the difference £ 2 ( s o) = z — z for fixed z . Notice that £ 2 ( s o) = (z^ — z^)(s) + (z^(s) - z( 2 \s)). 

Observe that the difference 8 = z^ — z^ satisfies the differential equation 

dS 



ds 



= Az« + O a . w (ez^) + SizW (e 2 ) + O s , w (\\z^\\ 2 ) = 



w. 



Here z' 1 ', and hence w, is to be regarded as a predetermined function for which we already have good 
bounds. Indeed, let § = {(so,s) : ^ sq < 2, sq ^ s ^ K} and || ■ stand for the C 3 norm on this set. 
According to (A. 2), \\w — Az^'Hc 3 ^ Ce 2 whereas, recalling that all eigenvalues of A are proportional to 
Ai, || A^r^J 11^3 < C|Ai|e. In other words, ||u>||c 3 . ^ G|Ai|e. As S(s ) = 0, we have 



(i) 



,(2) 



)(«) = 8(a) = f w(r)dT. 



Because s is C 3 in sq, it follows that IK-z^ 1 * 1 — z^)(s)||c3 ^ G|Ai|e. By (3.13b), the remaining contribution 



,(2) 



reads 



,(2) 



,,(o) 



dr. 



2(1 



We have seen above that ||s — s||c 3 ^ Ce, which implies \\z^ 2 \s) — z^- 2 \s)\\c3 ^ Ce 2 and finally ||^2( s o)l!c 3 ^ 
C|Ai|e. ■ 

Remark A. 4. It follows from the previous proof that, under the conditions of Proposition 3.1, 
(A.10) ||0, O S|| < C. 

Indeed, d Zo s = <9 Zo £i, as d Zo s = 0. From the fixed point equation £j = G(£i) we get d ZQ £,\ = (1 — 
DG(^i))^ 1 (<9 Zo G)(£i) and then the claimed bound. Moreover, 

(A.ll) ||a zo z-l|| <C|Ax|. 

For let z(s) = z(s; so, Zo) be the solution to (3.8b) with z(sq; sq, Zq) = Zq and recall that s depends on 
(s , z ). By definition, z — z(s; s , z ) so that d Zo z — d s z(s)d Za s + d Zg z(s) = 1 + 0(Ai) by the bounds in 
Lemma A. 3. 

Let us view the solution 

(A. 12) z(s) = z(s, s, z), z(s) = z 

to equation (3.7b) as a function of three variables and abbreviate d s = d/d s , di = djdzi, di 1 ...i k — di 1 ■ ■ ■ di k , 
and ds — d/ds. 

Proposition A. 5. Assuming ||z||/|Ai| is small enough, we have, for 0^fc + Z + m^3 and s ^ s, the 
following bounds: 

||<9™<^z( s )|| sC C\\z\\e^^ s - $s > 
dTd l A... lk z(s)\\ < ^L-e^t-') (k > 0). 

Proof. The initial conditions (dz / dz)(s) = 1, dijz(s) = 0, and dijkz(s) — follow from (A. 12), as the 
z-derivatives can be computed after evaluating z at s = s. Similarly, taking s-derivatives of 



s(s) = e^ A [i + J' e-^-^h^T, z(r)) dr^j 



yields first, analogously to how the identities below (A. 5) were obtained, 
dsz(s) = —Az — hi(s, z) 

d 2 s z(s) = A 2 z + A/ii(s, z) - d s hi(s, z) - Dh^s, z)d s z{s) 

9|z(s) = -A 3 z - A 2 h 1 (s, z) + 2Ad s h 1 (s, z) - d^h^s, z) + ADh^s, z)d s z{s) 

- D{d s h x ){s, z)d s z(s) - Dhx{§, z)^d s z(s) - D 2 hi{§, z){d s z(s), d s z(§)) 

as 

- Dh 1 (s,z)d 2 s z{s). 

These formulas can then be differentiated with respect to z in order to find higher-order initial conditions. 
As hi(s, z) = C(||z|| 2 ), we obtain the following estimates: 

\\d l s z(s)\\ ^ C\\z\\ 1 = 1,2,3 

\\d l s diz(s)\\ ^ C 1 = 1,2 

\\d s d ijZ (s)\\ <c. 

Combining (3.7b) and (A.l), 

d . z • Az + z ■ hi(s, z) , . ,,, , ... Ai , , 
^IN = ^ ~ ' ^ Ai||z|| + \\h x {a,z)\\ < y||z|| 

if ||z||/|Ai| is small enough. Below, ||z||/|Ai| will always be assumed small enough. Thus, for all s > s, 
(A.13) \\z{s)\\ ^ \\z\\e^( s - s \ 
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Differentiating (3.7b) with respect to various components of z, we obtain the variational equations 



(A.14) 


'i-diZ 
as 


= (A- 




z))diZ 




(A.15) 


d a 

dS d - z 


= (A- 




z))dijZ 4 


D 2 h 1 (s,z)(d l z,d 1 z) 


(A.16) 




= (A- 




z))d ijk z 


+ D 2 h 1 (s,z)(d l z 1 d jk z) 



+ D 2 h 1 (s, z)(d k z, dijz) + D 2 h 1 (s, z)(d j z 1 d ki z) + D 3 h 1 (s, z)(d l z,d k z,d ij z). 
Combining (A.f4), (A.f), and (A.f3), we have 

(A.17) \\diz(s)\\ ^ e^ a ~ s) 

in analogy with (A.13). Combining (A.15), (A.l), (A.13), and (A.17), 

Applying Corollary A. 2, 

(A.18) ||4-z( S )|| < ^^ (S ' S) - 

Similarly, combining (A.16), (A.l), (A.13), (A.17), and (A.18), 



C 



(A.19) \\di jk z(s)\\ < — 

Differentiating equations (3.7b), (A.14), and (A.15) with respect to s produces equations for d\z, d\diZ, 
and dgdijZ. For example, 

~d s z = (A + Dh 1 (s, z)) d$z 
as 

^-djz = {k + Dh 1 {s,z))d 2 z + D 2 h 1 (s,z){dsZ,d~ s z). 
as 

Such equations can be handled in a similar fashion and it is easy to verify that the additional s-derivatives 
do not change the bounds by more than a constant prefactor. The bounds with m = in the proposition 
are now clear. The bounds with m ^ follow immediately from the appropriate differential equation; for 
instance, bounding the right-hand side of (A.14) yields the bound on ^diZ. ■ 

Proof of Proposition 3.2. We first view the error terms £, k with 1 k 3 as smooth functions of (s, z) and 
bound their derivatives with respect to these variables. The bounds on the C 3 -norms with respect to so 
follow by the chain rule. Bounding the C 3 -norms of £4 and £5 is trivial and is done at the end of the proof. 
Throughout the proof, || • || C 3 will stand for the C 3 -norm with respect to sq. 

Terms £1 and £2. It is convenient to express £1 and £2 in the form 

/•OO fT 

£1=/ *>i(s + t) t P(s + t) / e^-^hxis + ^zis + ^dtdr 
Jo Jo 

/>oo 

£2 = / h 2 {s + t, z(s + t)) dr. 
Jo 

First of all, \\(b 1 P)(s + t)||c 3 ^C||&i~P||c 3 for every r, because p — sqWc 3 ^ C, so that 

/•OO PT 

<C||67P|| C 3 / / e x ^-^\\h 1 {s + i,z{s + i))\\ c ,didT 
Jo Jo 

||£ 2 ||c^ / \\h 2 (s + t, z(s + t))\\ c s dr. 
Jo 

Since hi (s,z) and h,2(s,z) are periodic in the variable s, their partial derivatives of any order (less than 
four) with respect to s are periodic functions of s and can be bounded exactly as /ii(s, z) and /i2(s, z). To 
save a considerable amount of space, we write 77 = s + £ and £ = (s + £, z(s + £)) below. 
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Notice that the first three total s-derivatives of z(rf) are 



— z(rj) = d s z(f])+d s z(r)) 
^z( V ) = d 2 s z( V ) + 2d s d s z(v) + d 2 s z(rj) 
^z{r,) = d 3 s z(v) + 3d 2 s d s z(v) + 3d s djz( v ) + djz{r,). 



Taking i-derivatives of the first two formulas above, 



-^diz(rj) = d s diz(rj) + d$diz(r)) 
as 

^ij z ( r l) = d s dijz(r]) + d$dijz(r)) 



ds 2 



diz{r)) = d 2 diz(i]) + d s d s diz(rj) + djdiz(rj). 



Proposition A. 5 then implies the bounds for k + I ^ 3: 



ds 1 



z(rj) 



< c \\z\\ e 



2 s 



ds 1 



c 



e 2 



H 



(k > 0). 



These will be used to bound the C 3 -norm of /ii(C). To this end, we compute 



ds 
d 2 



h 1 (Q=d s h 1 (C)+Dh l (<:)~z(r 1 ) 

ds 

= 0(\\z\\ 2 e x ^) 



ds 2 
ds 3 



ds ds z \ ds ds 



0(\\z\\ 2 e^) 



fc!(C) - 0?fci(C) + D{dlh x ){Q^z{r,) + 2D(d 2 s h 1 )(C)^z(v) + 3D(0 a fci)(C)-^z(»j) 

ds ds ds z 

+ 3D 2 (9 5 /m)(C) (jUfa), ±z(n)\ + DhxiQ^zfr) + 3D 2 ftx(C) (^Ufa). ^z{r,) 



d 



d 



+ D*h 1 (C)[~z( V ) ) ~z(r,),~z(rj) 



ds ds 



ds 



0(||z|| 2 e Al «) 



2.'i 



dzi 

d 2 

dzidzj 



d mo = ^MO^fa) 



= 0(||z||e^) 
MC) = Dh^Qdijziri) + D 2 h 1 (C)(d i z(r 1 ),d j z(r 1 )) 



"" h 1 (C)=Dh 1 (C)d ijk z( V )+D 2 h 1 (C)(d i z(r ] ),d jk z(r])) 



dzidzjdik 

+ D 2 h\{Q)(djz(rj),dkiz(rj)) + D 3 hi(C)0iz(r]),dkz(r]),d l: jz(ri)). 
Taking z-derivatives of ^fti(C)> ^2^i(C)i an d the resulting expression for ^z ds 

hi(C), we get 

''' MC) = Did.hjiQdiziri) + DhiiQ^-Aziv) + # 2 MC) ( £z(v), 



'/:,(/•- ,/ s V ds 

= 0(\\z\\e Xl6 ) 

''' r /i 1 (0- J D(a s 2 /ii)(C)^z(^) + 2i?(9 s /M)(C)^9^(r/)+2^ 2 (a^ 1 )(C) I ^z( V ), d iZ {r,) 



dzids 2 s ds V c£s 

+ 2D 2 mo + ^MC) jUfa). jUfa) 

= 0(||z||e^) 

<2 3 



diidz 



— h 1 (C)^D(d s h 1 )(C)d ll z( v ) + D 2 (d s h 1 )(C) (diz{ V ),djz(ri)) + Dh 1 {Q-d ij z{r ) ) 



ds w " J w v ^ V rf « 
+ £> 2 MC) + r> 3 MO ^z(»/). ( -^c,/).^z(,/.) 5 



We bound the derivatives of fi2(s + t, z(s + r)) in exactly the same way. 

Term £3. Setting u(t) = bi(r) T P(r) — E, we have £3 = 2- tt(T)e( r ~ s ' A dr. Using the facts that v is 
2-periodic, that its integral vanishes, and that A is negative definite, 



poo roo 00 / p2 \ 

J v(r)e {T -^ A dr = J v(s + T)e TA dT = J2l^J v(s + T)e TA dTj 

= (^j\(s + r)e TA d^j (l-e'T 1 

= (J* v(s + t) ( e ^ A - 1) dr) (1 - e 2 *)" 1 . 



e 2kA 



Hence, 



d k 



J™ v^e^-^ dr = ^ »<*>(J + t) (e TA - l) dr^j (l - e 2A ) _1 . 



Recalling that A is diagonal, we obtain for each value of k the upper bound 



(A.20) 



d k 
ds k 



(u(r) e ( T -^ A ) 



dr 



,(*0 



rA, 



1 dr (1 



„2A, 



,(fc) 



Incorporating (s , = 0) h-> (s, z). We set z = and denote (s, z) = Hk(so, 0). 

d fe Z 



As z = z(s) = z(s(s ), s , 0), we have £# = J^z(s(s ), s , 0). The bounds 



^Ce (1 < k 3) 



follow from the fact that s is a C 3 -function of Sq and the bounds in (A. 2). 
Since s and z are functions of s , for any function u = u(s, z), 



d 

-i—'< 
ds Q 

£_ 

ds 2 ' 

d 3 

ds 3 ' 



ds dzi 
-j—OsU + -7—o- Zi u 
dso dso 

d 2 s d 2 z, 
-j-2°su + -tyO^u + 
ds^ ds^ 



ds \ 



ds 0/ / 



, ds dz t dzi dzj 

) 9||it + 2— — % 4 u+ — -j^di^.u 
dso dso aso dso 



d 3 s n d 3 Zi „ „ ds d 2 s „ 
-r^dgu + -rg-o^u + 2 — — ^dssti + 2 
asg asg aso asp 



d 2 s d£ ds d 2 z.; 



ds 2 dso dso ds 2 , 



d 2 s d d 2 £i d 

dsg dso s ds 2 , dso z * 



ds 

dsn 



ds 



d 2 ^ dzj 
asp dsn. 

di, dz,- d 



2 — ^—d si .u+^^—dz.z.u. 
dso dso dsg * dso dso dso " l 3 



Here summation over repeated indices is understood and we leave it to the reader to expand the remaining 
sn-derivatives on the last line. Using the bounds derived earlier, we then get 

W^iOW <C\\zfe^ 

d 



ds 
d^ 
ds 2 
f_ 



MO 
hi(C) 



^C(||z|| 2 + £ ||z||)e Al « 
^C(||z|| 2 + £ ||z||+ e 2 )e Al « 
C | llzir -t- :\\z\ 



Similar bounds are obtained for h 2 (C)- We conclude that 

||8i||ca < C\\b[P\\c* 



e 2 e 2 



|Ai| 2 



|Ai| 2 



-2 || C 3 



^ C-, 



|Ai| 

||£ 3 || C 3 s? Ceo. 

The final inequality involving £i holds because is independent of a. 
Terms £4 and £5. Writing £4 in the form 

£4 = (s — s) / &o((l — T )s + ts) dr 
Jo 

and recalling that s and s are both C 3 functions of sn allows us to estimate 

||£ 4 ||c3 <C||s-s|| C 3 ^Ce. 
Proposition 3.1 was used here. Finally, by the same proposition, 

IIMc = m 2 (so,o),d)\\ c3 <*Ce||E||, 
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which finishes the proof. 



Lemma A. 6. For all sq and a, we have \\d Zo Soo(so, 0, <x) || > 0. 
Proof. Differentiating both sides of (3.11) with respect to Zo, we get 

(A.21) = b (s oo )d S0 s oo + (llQ e TA dr^j ,d Zo zJ+R, 

where 



R = -b (§)d z J + / (^°°( b i" p - + T V *r\ , d Zo z\ 

/ poo \ \ 2 

^ (bJPYis + T^dr) (d za s)j +J^d Z0 Z 



k ■ 



Because (A. 20) holds for any periodic, zero-integral function, the two integrals appearing in R are 0(a) in 
the limit Ai — > 0. Terms d Zo s and d Zo z are bounded by (A. 10) and (A. 11), respectively. Estimating d Zo £i 
and d Z[) £,2, we conclude that 

\\R\\=0{a) + O f ' 



^From (A.21) and (A.ll), we have 



(A.22) 



|Ai| 2 



SA- 1 (l + 0(A 1 )) + 0(a) + 



ere 



as Ai — > 0. Since SA 1 = (— SjA, )iLi > if jxn is sufficiently small, then the first term on the right side 
of (A.22) dominates and thus ||9 Zo Soo|| > 0. 
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